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ABSTRACT 

We present a spectroscopic study of Leo V, a recently discovered satellite of the Milky Way (MW). 
From stellar spectra obtained with the MMT/Hectochelle spectrograph we identify seven likely mem- 
bers of Leo V. Five cluster near the Leo V center (i? < 3') and have velocity dispersion 2At^\ km 
s~^. The other two likely members lie near each other but far from the center (i? ~ 13' ~ 700 pc) 



and inflate the global velocity dispersion to 3.7 



+i;4 km s" 



Assuming the five central members are 



bound, we obtain a dynamical mass of M = S.Stl l ^ lO^^o {M/Lv = 75tlf[M/Lv]Q)- From the 
stacked spectrum of the five central members we estimate a mean metallicity of [Fe/II]= —2.0 ± 0.2 
dex. Thus with respect to dwarf spheroidals of similar luminosity, Leo V is slightly less massive and 
slightly more metal-rich. Since we resolve the central velocity dispersion only marginally, we do not 
rule out the possibility that Leo V is a diffuse star cluster devoid of dark matter. The wide separation 
of its two outer members implies Leo V is losing mass; however, its large distance {D ^ 180 kpc) is 
difficult to reconcile with MW tidal stripping unless the orbit is very radial. 

Subject headings: galaxies: dwarf — galaxies: kinematics and dynamics — (galaxies:) Local Group 
— (cosmology:) dark matter — 



1. INTRODUCTION 

The recent discoveries of thir teen ultra-faint Milky 
Way (MW) satellites (e.g., [WiUman et aP 120051 : 
IZucker et aT1l2006l : iBelokurov et al.ll2007| ) have reshaped 
the census of Local Group galaxies and increased our 
ability to test cosmological models using objects avail- 
able in our own neighborhood. The new satellites, 
discovered primarily wi th data from the Sloan Digi- 
tal Sky Survey (SDSS; lYork et all [2OOOI ). extend the 
galaxy lumin osity function by three orders of magnitude, 
to ~ lO'^ LrT^ (iKoposov e tTal. 20 081). Spectroscop ic sur- 
veys Ce.g.. lMartin'eraLi2007, . .SinK)n fc Gehall2007l SG07 
hereafter) of their few red giants reveal that these sys- 
tems extend scaling relationships according to which the 
least luminous dwarf galaxies are the mo st dominated 
by dark matter a nd most metal-poor fe.g.. lMateolll998t 
iKirbv et al.l[2008h . 

Leo V, the most recent addition to the ensemble of MW 
satellites, was originally detected as a modest (~ 4g 
overdensity of stars in SDSS data (jBelokurov et al.ll200 
"Paper I" hereafter). Slightly deeper photometry from 
the Isaac Newton Telescope (INT) reveals a blue hori- 
zontal branch (BHB) and red clump in addition to the 
red giant branch (RGB; Figure 3 of Paper I). The appar- 
ent magnitude of its BHB implies Leo V has My ^ —4.3 
and lies at a distance of ^ 180 kpc, placing it at a (de- 
projected) dist ance of just ~ 20 kpc f rom another SDSS 
dSph, Leo IV (jBelokurov et al.ll2007l ). As the systemic 
line-of-sight velocities of Leos IV and V differ by just 
~ 40 km s^^, the two objects may represent the smallest 

•OBSERVATIONS REPORTED HERE WERE OBTAINED 
AT THE MMT OBSERVATORY, A FACILITY OPERATED 
JOINTLY BY THE UNIVERSITY OF ARIZONA AND THE 

SMITHSONIAN INSTITUTION. 

Electronic address: walker IQ) ast . c am . ac. ukl 

^ Institute of Astronomy, University of Cambridge, UK 

^ Department of Astronomy, University of Michigan, Ann Arbor 

^ Steward Observatory, The University of Arizona, Tucson, AZ 



companions among the known MW satellites. 

With this letter we present an initial spectroscopic 
study of individual stars in Leo V. Despite heavy contam- 
ination by the Galactic foreground, we are able to iden- 
tify seven likely members of Leo V. We use the stacked 
spectrum of these members to estimate Leo V's mean 
metallicity, and we use their spatial and velocity distri- 
butions to investigate Leo V's dynamical state and plau- 
sible dark matter content. 

2. OBSERVATIONS, DATA & MEMBERSHIP 

On 2008 May 28 and 30 we obtained medium- 
resolution spectra of 158 Leo V targets using the Hec- 
tochelle spectrograph at the MMT Observatory. We 
maximized the chance of observing Leo V members by 
choosing targets based on proximity to the center of Leo 
V and to the locus of its RGB (Figure 3 of Paper I). 
We extracted and calibrat ed spectra fol l owing the proce- 
dure described in detail by lMateo et al.l (|2008l ). The Hec- 
tochelle spectra span the wavelength range 5150 — 5300 
A, in which the most prominent absorption feature is 
the Mg-I/Mg-b triplet (MgT). We measure the fine-of- 
sight velocity of each star by cross-correlating its spec- 
trum against co-added Hectochelle spectra of radial ve- 
locity standards. We also measure a composite mag- 
nesium index, SMg, from the weighted sum of pseudo- 
equivalent widths f or each individual line of the MgT (see 
iWalker et al1[200l . 

We follow SG07 in modelling the velocity error as the 
sum of random and systematic components: the error in 

We deter- 



the 



■th 



measurement is ct; 



ran,« ' ^sys* 

mine the random errors, which correlate with S/N, using 
a bootstrap method. For the i*'' spectrum, we gener- 
ate 1000 artificial spectra by adding randomly generated 
Poisson deviates to the counts at each pixel. After mea- 
suring the velocity for each artificial spectrum, we equate 
'^ran i with the variance of the artificial distribution. We 
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then use repeat Hectochelle observations of 98 stars (in- 
cluding 11 Leo V targets) to determine the systematic er- 
ror common to all measurements. The value tJsys = 0.35 
km s~^ best satisfies our requirement that the distribu- 



tion of (Pi 



'ran,2 + 2cr.L) ^'"^ rescmblc a 



sys 



Gaussian with unit variance. By the same procedure, 
we estimate the systematic error of SMg to be 0.091 A. 
Among likely Leo V members in our sample, the mean 
(median) velocity error is ±2.2 km s^^ (±1.6 km s~^). 

Table [T] presents spectroscopic data for the 11 stars 
that have velocity in the range 162 < V/ (km s~^) < 
180, around the systemic mean reported in Paper I. The 
complete data set is available as an electronic table. 

Leo V's low luminosity implies that most of the stars 
overlapping its red giant branch are late-type dwarfs in 
the Milky Way foreground. To the eye, the Leo V popula- 
tion appears in Figure [T] as a clump with narrow velocity 
distribution at small radius and weak EMg. We quan- 
tify this separation using an expectation-maximization 
(EM) algorithm similar to those de vised for determ ining 
membership i n open clusters (e.g ., lSanderg|ll971h . Our 
algorithm fsee lWalker et al]|2009l for details) iteratively 
evaluates V and EMg distributions and assigns to each 
star a probability, Pm^ of belonging to the member pop- 
ulation. The algorithm converges after identifying seven 
stars as likely members of Leo V. 

Three stars (L5-31, L5-60 and L5-116 in Table [J) that 
would have been considered members after applying a 
conventional Bcr velocity threshold are identified as fore- 
ground {Pm < 0.01), based on strong EMg and large 
distance from the Leo V center. Two other stars (L5-52 
and L5-57) at large radius (i? ^ 13') are likely mem- 
bers (Pm = 0.7 ± 0.3 and Pm = 0.8 ± 0.3, respectively, 
where the reported error is derived from Monte Carlo 
tests in which we repeat the EM algorithm after resam- 
pling V and SMg from the sets of artificial spectra) based 
on their weak EMg values, which fall almost exactly on 
the mean value for Leo V. We find no members between 
3' < i? < 13'. Given Leo V's measured half-light ra- 
dius of Hi = 0.8' ± 0.1' (Paper I), the presence of two 
members beyond this gap is unexpected unless Leo V's 
morphology is highly distorted (see Section [5]). 

3. VELOCITY DISPERSION AND DYNAMICAL MASS 

Assuming Leo V has a Gaussian velocity distribution 
with mean (V) and variance dy , the data have 2D 

likelihood L{{V),<Jv,) « n^i[(< +^Fj"'/'exp[(y, - 
{V))'^/{2[ay +cryj)]*^' , where exponent M = 1 for mem- 
ber stars and M = for contaminants. The values Mi 
are unknown, but we can use the membership probabil- 
ities as weights in evaluating the expected log-likelihood 
given by^ 

N 



i?(lni) = -i^PM.ln«+4) 



1 ^ 

2 ^ 

i=l 



P 



Mi 



'Vo 



const. (1) 



* Our results are insensitive to our implicit assumption that 
the velocity error distributions are Gaussian. If we instead use 
the exact error distributions we obtain from the sets of artificial 
spectra, we obtain the same constraints on the velocity dispersion. 
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Fig. 1. — MMT/Hectochcllc spccotroscopic data for red giant can- 
didates in Leo V. Top: Sky positions of measured stars. Solid circles 
(crosses) indicate likely members (nonmembers) , for which the EM al- 
gorithm returns Pm > 0.5 {Pm < 0.5). An arrow points toward the Leo 
IV satellite, which has central coordinates (27, —165). Middle: Pro- 
jected distance from the Leo V center vs. line-of-sight velocity. Bottom: 
Magnesium index vs. velocity. The two members at large radius fall in 
the middle of Leo V's SMg distribution. 

Because we are interested primarily in ay^ , we follow 
iKleyna et alJ (j20Q4[ ) in considering the one-dimensional 
likelihood given by Lid«) - L[{V) , al^)d{V) . 

Again treating the probabilities Pm as weights, we eval- 
uate the ID expected likelihood 

E{\Ti Lid) = - \\na-]-{c- — 
I Z \ a 

1 ^ 

" 9 H ^A^' + O + const, (2) 



where a = E^i^M./(CTy. + cr^J, b = 

Ell PM^V./ia'y^ + <), and c = Eti huV.'/i^v, + 
ctvq)- We obtain error bounds by evaluating the area 
beneath the curve given by Lid — exp[_E(ln Li^)]. 
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TABLE 1 

Hectochelle Spectroscopy of Leo V (abridged: complete table is available in electronic version) 



Target 



L5-U01 
L5-002 
L5-004 
L5-007 
L5-008 
L5-031 
L5-052 

L5-055 
L5-057 
L5-060 
L5-116 



"2000 
(hh:mm:ss) 



<52000 
(dd:mm:ss) 



R 

(arcmin) 



(mag) 



g - r 
(mag) 



(km s"-*-) 



SMg 

(A) 



Pm 



11:31:13.21 
11:31:10.59 
11:31:02.88 
11:31:01.66 
11:31:15.91 
11:31:46.64 
11:31:14.11 

11:30:14.90 
11:31:04.15 
11:30:53.03 
11:32:22.01 



+02:12:51.6 
+02:14:09.5 
+02:13:14.7 
+02:11:25.2 
+02:10:57.6 
+02:09:10.9 
+02:26:36.4 

+02:14:30.1 
+02:26:56.5 
+01:58:20.2 
+01:56:49.3 



1.0 
1.7 
2.7 
2.7 
10.1 
13.4 



20.43 
19.85 
19.53 
19.67 
19.53 
20.80 
20.90 



0.78 
0.71 
0.74 
0.73 
0.70 
0.50 

0.59 
0.42 
0.77 
1.00 



174.8 
173.2 
168.8 
176.8 
175.2 
165.6 
166.0 
162.8 
179.2 
176.7 
174.6 



±0.9 
± 1.5 
± 1.6 
±2.1 
± 1.8 
±2.4 
± 2.6 
± 1.6 
±3.7 
±3.2 
±2.0 



0.07 ±0.38 
0.68 ± 0.16 
0.55 ± 0.16 
0.24 ±0.32 
0.27 ±0.22 
1.52 ± 0.19 
0.39 ± 0.25 
0.23 ± 0.20 
2.46 ± 0.32 
0.37 ±0.21 
1.67 ±0.52 
0.90 ±0.29 



"TUDTT 
1.000 
0.996 
0.994 
0.996 
0.000 
0.735 



± 0.000 
± 0.000 
± 0.003 
± 0.005 
± 0.002 
± 0.000 
± 0.353 



13.7 20.55 

13.8 21.31 
15.4 19.94 
24.4 21.00 



0.000 ± 0.000 
0.788 ± 0.304 
0.000 ±0.290 
0.000 ± 0.000 
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Fig. 2. — Mean velocity and velocity dispersion of Leo V, from 
the entire kinematic sample (~ 7 members; left panels) and from stars 
within 5' of the Leo V center (~ 5 members; right panels). Contours in 
upper panels enclose 68%, 95% and 99% of the volume underneath the 
likelihood surface described by Equation^ Lower panels give the ID 
likelihood of the velocity dispersion, with the shaded region marking 
the central 68% of the area under the curve. 

Because of the bimodal distribution of R among the 
Hkely Leo V members, we measure two sets of distri- 
bution parameters. Hereafter, "global" refers to the 
entire sample (7 members), while "central" pertains to 
a sample restricted to stars with i? < 5' (5 members). 

Figured (top) displays contours in the {{V),avo) plane 
that enclose 68%, 95%, and 99% of the volume under- 
neath the surface given by L = exp[£'(lnL)]. While the 
mean velocity and global velocity dispersion are well con- 
strained, we resolve the central velocity dispersion only 
marginally. This result holds also when we consider the 
ID likelihood of the velocity dispersion (lower panels in 
Figure [2|) . We measure a global velocity dispersion of 

avo = 3-7^^'^|^2 3) ^'^'i ^ central velocity disper- 

sion of — 2.4^^'^|]'[2 4^ km s^^, where errors indicate 

68% (95%) confidence levels. Thus at 95% confidence 
we rule out zero dispersion only for the global sample 
and cannot state conclusively that we resolve the central 
velocity dispersion. 



Assuming the five central members are bound, we esti- 
mate the dynamical mass of Leo V using a crude "mass- 
foUows-light" (MFL) model in which the density profile 
of any dark halo is proportional to that of the stellar 
component, and which implies total mass M = rjriiay^ 
(jlllingworthl Il976f ). For the purpose of comparing to 
published estimates for other satellites we adopt 77 = 
850Mopc~^km~^s^, the value corresponding to a iKingI 
(|1962( ) profile with co ncentration p arameter characteris- 
tic of brighter dSphs (|Mateolll998f ). With this value and 
the measured central velocity dispersion we obtain M 



3.3 



+9.1(+46) 



75 



+230(+1200) 



2.5(-3.3) ^ lO^Afo and MmFI^/Lv - . -_58(_74) 

While the lower mass limit should be interpreted as con- 
sistent with a purely stellar population free of dark mat- 
ter, even the upper limit gives Leo V one of the lowest 
dynamical masses of the SPSS satell i tes with measure d 
kinematics (SG07. lMartin et al.l[2007HGeha et al.ll2008^ . 

4. METALLICITY 

We est imate meta.llicity by comparing spectra to the 
library of' Lee et all (|2008h . based on a regularly sampled 
grid on which [a/Fe] = 0.4 covers a metallicity range 
-4.0 <[Fe/H]< -0.5, effective temperature 3500/^ < 
Tcff < 9750K, and surface gravity 0.0 < logg < 5.0. 
Since individual LeoV spectra at i? = 25000 have in- 
sufficient S/N to perform a reliable spectral analysis we 
averaged the 5 central member spectra and, after suitable 
Gaussian smoothing, rebinned this average spectrum to 
a resolution of i? = 10000. After continuum normalizing 
both data and model spectra, we compare the data di- 
rectly to a suitable subset of the synthetic spectra using 
a masked least-squares minimization. The adopted mask 
isolates regions where significant absorption lines (typi- 
cally with peaks < 0.97 of the continuum) appear in the 
model spectrum (Figure [3]). 

This free-form fit gives a well-defined solution centered 
on Tcff = 5000 ± 250K, [Fe/H]= -2.0 ± 0.3 dex and 
log 5 = 2.0 ± 1.0. Reassuringly, if we use the relatio nship 
logio(reff) = 3.877 - 0.26(.9 - r) (jlvezic et al.ll2006[ ). the 
average color {{g — r) — 0.72) of the five central members 
independently implies (Toff) = 4900K, consistent with 
the temperature we derive directly from the continuum- 
normalized spectra. Anchoring the value of Teff for the 
spectral fit tightens constraints on gravity and metallic- 
ity, giving logg = 2.0 ± 0.5 and [Fe/H] ^ -2.0 ± 0.2 dex. 
The average color, magnitude, surface gravity and effec- 
tive temperature are fully consiste nt with earl y K-giants 
at the distance of Leo V (see, e.g.. lGravll2005l . p. 57). 
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Fig. 3. — Top: Stacked spectrum (continuum- normalised, aver- 
aged and rebinned) for the central 5 Leo V members. Bottom: the 
best fit model spectrum; and in blue above the mask used in the 
minimisation. 

5. DISCUSSION 

Because we do not necessarily resolve the central veloc- 
ity dispersion, we can neither place a meaningful lower 
limit on the dynamical mass nor state definitively that 
Leo V is a dark-matter-dominated dSph rather than a 
diffuse star cluster devoid of dark matter. In fact the 
data lend some support to the latter interpretation, sug- 
gesting that Leo V is an outlier with respect to two em- 
pirical scaling relations among dSphs. First, its dynami- 
cal mass-to-light ratio (M/Ly = TStgf [M/Ly]©) likely 
falls short of the more extreme values (~ 10'^~^[M/Lv]q) 
exhibited by dSphs of similar luminosity (e.g., Figure 15 
of SG07). Second, at [Fe/H]= -2.0 ± 0.2 dex, Leo V 
is metal-rich compared to dSphs of similar luminosity, 

which typically have [Fe/H] 2.5 (SG07, iKirbv et al.l 

[2008h . 

Despite ambiguities regarding the classification of Leo 
V, one conclusion is clear: the available data are incon- 
sistent with the notion that Leo V is an equilibrium sys- 
tem with half-light radius rh = 0.8' ± 0.1' (~ 50 pc), 
the value measured from the surface brightness of its red 
giants (Paper I). Assuming Leo V is a uniformly sam- 
pled (c.f. Figure [T]) Plummer sphere, the probability of 
finding at least two members, in a sample of seven, at 
R > lOrh is ~ 10"''. Thus our detection of two members 
at R ^ 13' (~ 700 pc) implies a distorted morphology of 
the sort that might arise if the two outer members be- 
long to an unbound tidal stream. The extreme velocities 



of these two stars — they contribute the smallest and the 
largest velocities among members — further suggests that 
they are unbound. 

Tidal interactions with the MW produce stellar 
streams observable in several MW satelli tes — for exam- 
ple, the Pal 5 globular cluster ( Odenkirchen et al.ll2003fl 
and t he Sagittarius dwarf galaxv (e.g.. iMaiewski et all 
|2003() . However, Leo V, at a distance of D 180 kpc 
(Paper I) , will be unaffected by MW tides unless its orbit 
is strongly radial. If Leo V is on a radial orbit, then its 
journey from near the center of the MW to its present 
position must have taken 1 Gyr. Assuming the two 
outer stars were once bound to Leo V and were stripped 
at the most recent pericentric passage, then they could 
have plausibly moved a distance of ^ 700 pc along the 
tidal tail in this time (see Equation 7 of Odenkirchen et 
al. 2003). 

We argue in Paper I that Leos IV and V are likely 
companions, in which case their systemic velocities im- 
ply a nearly circular orbit that would render MW tides 
unimportant. In this case the two satellites might inter- 
act with each other in filamentary substructure of the 
sort produced in cos mological N-body simulations (e.g., 
iDiemand et al.ll2007l ). The fact that Leo V's two outer 
stars lie along the line connecting Leo IV to Leo V (Fig- 
ure[Tl top) provides some support for this scenario, as do 
the extended distributions of BHBs around both systems 
(Fig. 4 of Paper I). 

Alternatively, Leo V may be a loosely bound star 
cluster losing stars to evaporation. It is not clear 
whether such a large cluster could form in isolation, 
but Leo V may have been stripped from a progenitor 
among the more luminous MW dSphs. The best can- 
didate is Leo II {D - 220 kpc, VbsR 22 km s^^), 
which lies within ~ 1.1° of the orbit of the prospective 
Leo IV/V stream. It is also possible that a progeni- 
tor with low surface brightness lurks outside the SDSS 
footprint, in which case it should be detectable with 
data from upcoming deep-in iaging surveys such as Pan- 
STARRS (Kaiser et al."2002^ and the Southern Sky Sur- 
vey (Keller ct al. 2007). 
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